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Abstract: We report a mathematically rigorous technique which facilitates 
the optimization of various optical properties of electromagnetic fields. The 
^^ ■ technique exploits the linearity of electromagnetic fields along with the 

CS| \ quadratic nature of their interaction with matter. In this manner we may 

decompose the respective fields into optical quadratic measure eigenmodes 
(QME). Key applications include the optimization of the size of a focused 
spot, the transmission through photonic devices, and the structured illu- 
mination of photonic and plasmonic structures. We verify the validity of 
the QME approach through a particular experimental realization where the 
size of a focused optical field is minimized using a superposition of Bessel 
beams. 
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1. Introduction 

The decomposition of fields into eigenmodes is a well established technique to solve various 
problems within physical sciences. The most prominent example refers to Schrodinger's equa- 
tion, within the field of quantum mechanics, where energy spectra of atoms are determined via 
the eigenvalue spectra and associate wavefunctions of the Hamiltonian operator. Indeed, elec- 
tron orbits are eigenmodes of the energy, angular momentum, and spin operators CI and as 
such they deliver fundamental insights into the physics of atoms. Within classical mechanics, 
modes of vibration of music instruments give, for example, their resonant frequencies while 
their spectrum is associated with the shape of the instrument |2|. In the optical domain, mode 
decomposition is used in order to describe light propagation within waveguides 13], photonic 
crystals [4], and optical cavities lEl. In the case of waveguides and photonic crystals, for ex- 
ample, eigenmodes describe electromagnetic fields that are invariant in their intensity profile as 
they propagate along the fibre or crystal. Additionally, these modes are orthogonal and as such 
light coupled to one of these modes remains, in theory, in this mode forever. This optical mode 
decomposition can be expanded to include additional operators such as orbital and spin angular 
momentum |6|. 

In this paper, we report a novel method which we term "Quadratic Measure Eigenmodes 
(QME)" which represents a generalization of the powerful concept of eigenmode decomposi- 
tion within the field of optics. Crucially, it is shown that eigenmode decomposition is applicable 
to the case of any quadratic measure which is defined as a function of the electromagnetic field. 
Prominent examples of optical quadratic measures include the energy density and the energy 
flux of electromagnetic fields. The QME method makes it possible to describe an optical system 
and its response to incident electromagnetic fields as a simple mode coupling problem and to 
determine the optimal "excitation" for the given measure considered. Intuitively, a superposi- 
tion of initial fields is optimized in a manner that the minimum/maximum measure is achieved. 
For instance, the transmission through a pinhole is optimized by maximizing the energy flux 
through the pinhole. 

From a theoretical perspective, the QME optimization method is mathematically rigorous 
and may be distinguished from the multiple techniques currently employed ranging from ge- 
netic algorithms |7| and random search methods |8| to direct search methods 0. The major 
challenge encountered in any such approximate optimization and engineering of optical prop- 
erties is the fact that electromagnetic waves interfere. As such the interference pattern not only 
makes the search for an optimum beam problematic but crucially renders the superposition 
found unreliable, as the different algorithms may converge on different local minima which 
are unstable with respect to the different initial parameters in the problem. In contrast, our 
proposed QME method yields a unique solution to the problem and directly determines the 
optimum (maximal/minimal) measure possible. 

In the first part of the paper, we introduce the background of the QME theory and show its 
properties in a general context of optimizing the quadratic measures of interfering waves. In 
the second part, we apply the QME formalism to maximize the transmission through apertures 
and to minimize the focal spot size. For these applications, we describe, respectively, the elec- 
tromagnetic field as a superposition of scalar Laguerre-Gaussian beams and vectorial Bessel 
beams. In the final part of the paper, we report a particular experimental implementation of the 



QME method using computer controlled spatial light modulators to squeeze the spot size of a 
superposition of Bessel beams. The paper concludes with a discussion of the particular results 
obtained and with general comments on the versatility of the QME method to a wide range of 
problems. 

2. Fundamental concepts 

The QME method is based upon two fundamental properties of the electromagnetic field and 
its interactions. Firstly, the approach relies on the linearity of the electromagnetic fields, i.e., 
the sum of two solutions of Maxwell's equations is itself a solution of them. As we consider 
free space propagation, this criteria is satisfied. The second property relates to the interaction of 
the electromagnetic field with its environment. All such interactions can be written in the form 
of quadratic expressions with respect to the electric and magnetic fields. Examples include 
the energy density, the energy flow, and Maxwell's stress tensor. This allows us to designate 
appropriate QME to various parameters (e.g. spot size) and subsequently ascertain the optimal 
eigenvalue which, in the case of a spot size operator, yields a sub-diffraction optical spot. In 
this section, we present the details of the theory underpinning our approach. 

2.1. Theoretical background 

Electromagnetic waves To demonstrate our method, we consider monochromatic solutions 
of the free space Maxwell's equations: 

V-eoE = 0, V-;UoH = 0, V xE= -ijioCoU, VxH = ieo(oE, (1) 

where E and H are the spatial part of the electric and magnetic vector fields and where Eq and /io 
denote the vacuum permittivity and permeability. The time dependent carrier wave is given by 
exp(i(jL)^). These monochromatic solutions of Maxwell's equations can be written in an integral 
form linking the electromagnetic fields on the surface A with the fields at any position r. 



^iir)= f Pijiry).^jir')dS' 

JA 



(2) 



where \/2^ = (y'eoE,y7^H) is a shorthand for the two electromagnetic fields having six 
scalar components ^t. The integration kernel Vij corresponds to a propagation operator giving 
rise to different vector diffraction integrals such as Huygens, Kirchhoff |10|, and Stratton- 
ChuEl. 

Quadratic measures Crucially all "linear" and measurable properties of the electromagnetic 
field can be expressed as quadratic forms of the local vector fields and are therefore termed 
quadratic measures. For instance, the time averaged energy density of the field is proportional 
to ^* • ^ = l/2(eoE* • E + ^qW • H) while the energy flux to 1/2(E* x H + E x H*). The 
asterisk * stands for the complex conjugate. Integrating the first quantity over a volume deter- 
mines the total electromagnetic energy in this volume, and integrating the normal energy flux 
across a surface yields the intensity of the light field incident on this surface. All the quadratic 
measures m^ can be represented in a compact way by considering the integral 

m^ = j ^^Kij^jdr = {^\ K\^)v (3) 

where the kernel Kfj = Kj- is Hermitian with t the adjoint operator including boundary effects 
for finite volumes. Table [T] enumerates some operators associated to common quadratic mea- 
sures. The integrand part of most of these quadratic measures corresponds to the conserving 



densities, which together with the associated currents are Lorentz invariant JS]. The volume, 
over which the integral is taken, does not need to be the whole space and can be a region of 
space, a surface, a curve, or simply multiple points. To account for this general integration 
volume, we broadly term it the region of interest (ROI) in the following. 



Operator 


^^t^ij^i 


EO 


eoE*-E + AtoH*-H 


10 


(E*xH + ExH*)-e^ 


SSO 


r2(E*xH + ExH*)-e^ 


LMO 


£oE*-(/5,)E + AioH*-(/5,)H 


OAMO 


eoE* • (/r X V)^E + MoH* • (/r x V)^H 


CSO 


i(E*-H-H*-E) 



Table 1. Common quadratic measure operators including the energy operator (EO), inten- 
sity operator (lO), spot size operator (SSO), linear momentum operator (LMO), orbital 
angular momentum operator (OAMO), and circular spin operator (CSO). The vector oper- 
ators include the subscript k indicating the different coordinates and e^^ the associated unit 
vectors. 



Quadratic measures eigenmodes Finally, using the general definition (O of the quadratic 
measure it is possible to define a Hilbert sub-space, over the solutions of Maxwell's equa- 
tions, with the energy operator defining the inner product. Furthermore, any general quadratic 
measure defined by (|3]) can be represented in this Hilbert space by means of its spectrum of 
eigenvalues and eigenfunctions defined by: 



A^Kr2) = / K-,-y(r2,ri)^XriMri = /c|^)y. 



Depending upon the kernel Kij or operator /c, the eigenvalues A form a continuous or discrete 
real valued spectrum which can be ordered. This gives direct access to the solution of Maxwell's 
equations with the largest or smallest measure. The eigenfunctions are orthogonal to each other 
ensuring simultaneous linearity in both field and measure. 

In the following, we study the case of different quadratic measure operators and their spec- 
tral decomposition into the QME. The convention for operator labeling we adopt is to use the 
shorthand QME followed by a colon and a shorthand of the operator name. In the following, 
we discuss some examples of different quadratic measure operators. 



QME: Intensity operator (QME:IO) The quadratic measure corresponding to the QME:IO 

measures the electromagnetic energy flow across a surface ROI: 
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\ f (E*xH + ExH*)-n(i5 
2 jRoi 



(4) 



where K is chosen such that it corresponds to the projection of the Poynting vector on the normal 
n to the surface. The eigenvector decomposition of this operator can be used to maximize the 
optical throughput through a pinhole or to minimize the intensity in dark spots. Considering a 
closed surface ROI surrounding an absorbing particle, the QME approach gives access to the 
field that either maximizes or minimizes the absorption of this particle. 



QME: Spot size operator (QMEiSSO) One way to define the spot size of a laser beam is 
by measuring the second order intensity moment (SOIM) w lfT2l . w can be expressed as 



(5) 




with m(^) as defined in Eq. (H) and the QMEiSSO defined by 

m^^^ = \ I |r-rop(E*xH + ExH*)-n(i5, 
2 Jroi 



(6) 



where r is the position vector and ro the centre of the beam. Accordingly, the square root of 
the QME:SSO eigenvalues multiplied by two gives direct access to the respective beam size 
provided that the intensity is normalized to one within the ROI, i.e., mS^' = 1; 



QME: Optical force operator (QMEiOFO) The optical force acting on a scattering object 
can be calculated by considering the momentum flux, given by Maxwell's stress tensor, across a 
surface ROI, surrounding the scattering particle. There are three quadratic measures associated 
with the optical force acting on the particle, one for each orthogonal direction: 
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(7) 



where m^^^ = (nix •,'^x , ^x ) corresponds to the measure of the force. The eigenvector 
decomposition of this operator can be used to maximize the optical scattering and optical trap- 
ping force on microparticles. This approach can be extended through the use of the angular 
momentum operator, defined locally by ir x V. 

3. Applications 

In this section, we put the QME concept into practice and provide a couple of examples which 
demonstrate the striking applicability of the QME formalism to answer different questions 
within the field of optics. One such question is determining the largest intensity that may be 
transmitted through a given optical structure such as metallic apertures 1 13|. Another question 
is what is the smallest optical spot one may achieve. Before we discuss these questions, we 
show how the operator formalism presented in section [2l is rendered into a practical formalism 
which we have also applied in the experiments described in section IH 

3.1. Practical implementation of the QME concept 

For practical purposes, in particular in terms of the experimental realization of the QME con- 
cept, the optimization procedure is spatially separated; that is we consider 1) an initial plane at 
the propagation distance z = Z\ where we can superpose a set of / = \...N fields E/(x, j,z) and 
H/(x, j,z) (N > \) shaped both in amplitude and phase and 2) a target plane at the propagation 
distance z = Z2 where we have the ROI within which the optimization is actually carried out. 
Due to linearity a superposition of fields in the initial plane is rendered into a superposition in 
the target plane both featuring the same set of superposition coefficients: 



E(x,j,zi) = ^a,-E,-(x,3;,zi) 

i=\ 

N 

H(x,j,zi) = ^(3,-H,-(x,3;,zi) 



i=\ 
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Based on this, the QME:IO defined in Q can be rewritten as 

mW=a*MWa. (9) 

The matrix M^^^ is slN x N matrix with the elements given by the overlap integrals 

mJ) = i / (E*(x,);,Z2) X Uj{x,y,Z2)^M^.y.Z2) x H*(x,3;,Z2)) -ndS. (10) 

This matrix is equivalent to the QMEiIO on the Hilbert subspace defined by the fields 
{E/(x,j,Z2),H;(x,j,Z2)}. M^^) is Hermitian and positive-definite which implies that its eigen- 
values Xj^'{k=\...N) are real and positive and the eigenvectors v^ are mutually orthogonal. 
Accordingly, the largest eigenvalue Amax = max iXj^ M and the associated eigenvector Vmax will 
deliver the initial plane (z = Zi) or target plane (z = Zi) superposition 

Emax (X, y,z) = Y. W,/E/ (X, J, z) , Hmax (x, J, z) = ^ %ax,/H/ (x,y,z)\ , (11) 

which maximizes the intensity within the ROI. We remark, that the case of a superposition of 
scalar fields Ui{x^y^z), the matrix operator ([TOl) takes on the simpler form 

^ij = / ul{x,y,Z2)uj{x,y,Z2)dS, (12) 



an approach used in section 1X21 to maximize the transmission through a pinhole. 
Similar to the QMEiIO, the QME:SSO defined in Eq. (O can be rewritten as 

m(^) = b*M(^)b, (13) 

where M^^^ and b must be represented in the intensity normalised base 

N yf] ^ N yf) \ 

^k{x,y,Z2) = Y, -^=-M^^y^Z2), Hk{x,y,Z2) = ^ -yd=-lii{x,y,Z2) > (14) 

^=1 y Af ^ ^=1 y Af ^ J 

in order to fulfill the requirement of unity intensity within the ROI (m^^^ = 1, see paragraph 
"QME: Spot size operator (QME:SSO)" in section [2Jl). M^^) is a A^ x A^ matrix with the ele- 
ments given by 

M^f = \ J^^^ |r-rop (e*(x,);,Z2) x Uj{x,y,Z2)^Mx^y.Z2) x n){x,y,Z2)) -ndS. (15) 



After transforming back to the original base we denote the eigenvalues of M^^^ as A^^ ^ and 

(2) (2) 

the eigenvectors as v^ \ Then, the eigenvector v^^^ associated with the smallest eigenvalue 
^min ~ ™^ ( \ ) corresponds to the smallest spot achievable within the ROI. The respective 
superposition is 

{ ^rmn{x,y,z) = £ v|^^^ .E,-(x,3;,z), Hrnin{x, y,z) = £ v|^^^ .H,-(x,3;,z) > , (16) 

[ i=l i=l J 

for z = Z\ and z = Z2- We employed this vectorial definition to minimize the size of a focused 
spot in section [34l using a superposition of vector Bessel beams. The spot size minimization 
performed both in section [33l using a superposition of LG beams and in the experimental sec- 
tion |4] is based on the simpler scalar version ot the QME: S SO matrix defined as 

Mif=[ \r-ro\^ir^{x,y,Z2)uj{x,y,Z2)dS. (17) 



3.2. Maximizing transmission through apertures using Laguerre Gaussian beams 
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Fig. 1. 2D intensity profiles together with the ROI highUghted by the red circle. The trans- 
mittance through the ROI is described by the coefficient T. (a) Transversal and (b) longitu- 
dinal cross section of a Gaussian beam (L = and P = 0). (c) and (d) QME superposition 
delivering the highest intensity in the ROI (R = wq). 



In this subsection, we consider a superposition of LG beams propagating in the z-direction 
and modulating the carrier wave Uc = exp(— i(A:oZ — cot)). In cylindrical coordinates, LG beams 
are defined by lfT4ll : 



Up{r,(^,z) = —j-r- r- / X 
q{z) \V2q{z) 



\L\ 



-q%z) 



r\L\ 



2r^ 
w\z) 



exp 



2q[z) 



iL0 



(18) 



with Zr = ^oWq/2, q{z) = z + kr, and w^{z) = 1 -\-z^/z^ where wq, ko, CO are the Gaussian 
beam waist, vacuum wave vector, and optical frequency, respectively. This beam profile is a 
solution of the paraxial equation for integer values of P and L corresponding respectively to the 
radial and azimuthal index of the beam. Within the paraxial approximation, the intensity of the 





Fig. 2. (a) Total transmittance through the ROI for the QME intensity optimized beam 
as a function of the ROI relative radius R/w{) and for different numbers N of LG modes 



considered, (b) Relative intensity 
QME intensity optimized beam. 
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of the LG modes (P = O..A^) decomposing the 



beam transmitted through a planar ROI defined by a disk centered on r = is proportional to 
/o^2;rr|w^(r,0,z)P<ir where 7? is the radius of the disk. The coefficient C^ = ^J P\ /[p ^ \L\)\ h 
the normalization factor such that the total intensity of the beam, for an infinite ROI, is unity. 

The transmission is maximized using the practical implementation of the QME concept de- 
scribed above in section [3Jl that is the QME:IO was assembled according to Eq. (O using 
the respective representation in cylindrical coordinates. We only considered the radial family 
of LG beams (L = 0) and performed the optimization in the plane at z = 0. Eigure[T] shows the 

final superposition u^^^{x^y^ 0) = Y!p=\ ^max p^pi^^y^ ^) i^ the case of a ROI with a radius equal 
to the waist of the Gaussian envelope. 

In Eig. [21 we observe that the maximal transmission achievable via the QME intensity opti- 
mized beams depends on the number of LG beams considered in the superposition and on the 
size of the ROI. Indeed, the QME for a smaller ROI needs a larger number of LG modes to 
achieve 100% transmission. 

3.3. Smallest focal spot using Laguerre Gaussian beams 

Using a superposition of LG beams we have also minimized the size of a focal spot using 
the representation of the QME:SSO ([Tt]) in cylindrical coordinates. It is important to note at 
this point that we only retain the intensity eigenmodes whose eigenvalues are within a chosen 
fraction of total intensity. This is equivalent to considering only beams that have a significant 
intensity contribution in the ROI. Intuitively, the optimization procedure may be performing 
so well that a spot of size zero is finally obtained if no intensity threshold is applied. EigureO 
shows the smallest spot superposition where we observe the appearance of sidebands just out- 
side the ROI. These sidebands are a secondary effect of squeezing the light below its diffraction 
limit. It is these sidebands that decrease the efficiency of the squeezed spot with respect to the 
maximal possible intensity in the ROI as calculated via the QME:IO. Using the ratio between 
these two intensities we can define the intensity Strehl ratio ifTSl for the QMEiSSO (see Eig. lib). 
We remark that both, the spot size and the Strehl ratio, show resonances as a function of the 
ROI size. This can be explained considering the number of intensity eigenmodes used for the 
spot size operator. Indeed, as the ROI size decreases, so does the number of significant inten- 
sity eigenmodes. Each time one of these modes disappears (step in Eig. |4]), we have a sudden 
increase in the minimum spot size achievable accompanied with an enhanced Strehl ratio as 
we drop the most intensity inefficient mode. Overall, the Strehl ratios determined in our studies 
predominantly exceeded values of 1 % even when spots were tightly squeezed. Therefore, the 
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Fig. 3. (a) Transversal and (b) longitudinal 2D intensity cross sections of the QME super- 
position delivering the smallest focal spot in the ROI (R = X) considering 25 LG modes. 
w/wo is the relative SOIM measured according to Eq. (|5]). The Strehl ratio in (a) is 4.5%. 
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Fig. 4. (a) Spot size as a function of the radius of the ROI for different number of LG 
modes considered. The right hand scale and filled curve indicate the numbers of intensity 
eigenmodes N^^' fulfilling the intensity criteria for the A/^ = 11 case. The arrows indicate 
the corresponding scales, (b) Ratio between the ROI intensity of the smallest spot size 
eigenmode and the largest intensity achievable in the ROI (Strehl ratio). 



observed decrease of intensity is not to severe in terms of potential applications of squeezed 
beams for optical manipulation and imaging. 

Analogously to the approach to calculate the smallest spot size in a planar ROI, we can 
determine the QME:SSO for a volume. In this case, the sidebands appear in the outside of the 
ROI in both the lateral and longitudinal directions. The different volumes considered in Fig. [5] 
suggest that there is an intrinsic link between squeezing light in the lateral direction and in the 
longitudinal direction. 

On a final note, we remark that squeezing light below its diffraction limit may be associated 
with the effect of super-oscillations |[T6l . This refers specifically to the ability to have a local 
/:- vector (gradient of the phase) larger than the spectral bandwidth of the original field. To visu- 
alize this effect, in the case of QME spot size optimized beams, we have calculated the spectral 
density of the radial wave- vector for the smallest planar spot fVl\. As shown in Fig.[6l this spec- 
tral density clearly identifies a spectral bandwidth (white background in Fig.[6l). Regions of the 
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Fig. 5. (a,c) Transversal and (b,d) longitudinal intensity cross sections of the QME:SSO for 
two different 3D ROL Strehl ratio: (a) 5.2%, (c) 1.5%. 



beam which exhibit locally larger wave- vectors than the ones supported by this spectral band 
width correspond to super-oscillating regions. The local wave vector is defined as 5^arg(w(r)) 
where arg(w) defines the phase of the analytical signal u. We observe, that super-oscillations 
occur only in the dark region of the beam. Additionally, when the ROI is large compared to the 
Gaussian beam waist wq, there are no super-oscillating regions. These only appear when the 
beam starts to be squeezed. 

3.4. Smallest focal spot using Bessel beams 

The paraxial approximation employed above in the case of LG beams can be used to describe 
sub-diffracting beams but breaks down when beams are tightly focused. As a consequence we 
must consider full vectorial solutions of Maxwell's equations. Here, we have chosen Bessel 
beams as a base-set and determined the superposition of Bessel beams which minimized the 
spot size in a planar finite ROI. Note that the problem of the finite intensity of Bessel beams ifTSll 
is easily circumvented here due to the finite ROI size considered. The monochromatic electric 
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Fig. 6. (a) Radial wavevector spectral density. Yellow highlights regions outside the spectral 
bandwidth, (b) Transversal cross section of the QME spot size optimized field intensity with 
yellow showing super-oscillating regions. 




Fig. 7. Intensity cross sections: (a) Airy disk for the maximum numerical aperture consid- 
ered NA= sm(Ofnax) = 0.1. The yellow dashed circle shows the position of the smallest 
zero-intensity circle taken as the ROI inside which the spot size is calculated. The spot size 
is normalized to the spot size of the reference Bessel beam, (b) Reference Bessel beam 
corresponding to the largest cone angle Omax- The SOIM of the reference Bessel beam is 
denoted as wb- (c) QME spot size optimized beam for a superposition of Bessel beams 
(G e [0, Omax]) for a large ROI highhghted by the dashed yellow circle. Strehl ratio: 2%. (d) 
QME spot size optimized beam for a small ROI. Strehl ratio: 0.2%. The gray-scaled region 
shows the sidebands while the color range the ROI. Notice that the two scales are different. 



vector field of the vectorial Bessel beam may explicitly be expressed as l|T9ll 



E 



Eocxp{iL(l)^iktz)) {aex^Pey)JL{ktr) 



ikt 



((a + /j3)exp(-/0)/L-i(fcr)-(a-/j3)exp(/0)/L+i(fcr))e^ 



(19) 



where kt =kosin{6) and k^ = ^ocos(O) are the transversal and longitudinal wave vectors with 
the characteristic cone angle of the Bessel beam, e^, e^; and e^ are the unit vectors in the 
Cartesian coordinate system. The parameter L corresponds to the azimuthal topological charge 
of the beam while a and j3 are associated with the polarization state of the beam. The magnetic 
field H was deduced according to Maxwell's equations and the QME: 10 and QME: S SO oper- 
ators are assembled according to the expressions (fTOl) and ([HI) , respectively. Figure [7] shows a 
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Fig. 8. (a) Relative spot size Ar/wB of the Bessel beam superposition as a function of the 
relative ROI radius R/Rb- The SOIM w^ and the ROI radius R^ are associated with the 
reference Bessel beam shown in Fig. [3b), where the ROI is indicated as dashed circle. For 
comparison, the red dot indicates the location of the reference beam in the Ar/w^ vs. R/Rb 
plot, (b) Strehl ratio vs relative ROI radius R/Rb • 



comparison between Airy disk, Bessel beam, and QME optimized spot considering a numerical 
aperture of NA= 0.1. As in the case of the LG beams, squeezing the focal spot is accompanied 
by side bands and a loss in efficiency shown by the Strehl ratio (see Fig. [8]). 

4. Experimental QME 

4.1. Experimental implementation of the QME concept 
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Fig. 9. Experimental setup. FP = focal plane, L = Lens. Focal widths: /i = 50 mm, 
/2 = 500 mm, h = fA = 400 mm, /s = 1 m. Laser: JDS Uniphase HeNe laser, Pmax = 
10 mW, A = 633 nm, SLM: Holoeye HEO 1080 P dual display system, resolution = 
1920 pixel X 1080 pixel, display size = 1 in x 0.7 in. CCD camera: Easier pilot piA640- 
210gm, resolution = 648 pixel x 488 pixel, pixel size = 7.4 p^m x 7.4 \xni. 



According to the theoretical foundations of the QME concept, the successful experimental 
implementation within the field of optics requires a linear optical system along with the ability 
to shape laser fields in both amplitude and phase. We have achieved this by using the exper- 
imental setup shown in Fig. [9l A HeNe laser beam is expanded and subsequently amplitude 
modulated by a spatial light modulator (SLM) display operating in conjunction with a pair of 
crossed polarizers. Analogously to a liquid crystal display on a computer or laptop monitor, the 



liquid crystal SLM display rotates the polarization of the incident light by an angle depending 
upon the voltage applied to the display pixels. The amplitude modulated beam is then imaged 
onto a second SLM display through a pair of lenses. This second SLM display along with a sub- 
sequent Fourier lens and aperture served to modulate the phase of the laser beam in the standard 
first order configuration ll20l . The field modulations of interest were encoded as RGB images 
where the blue channel represented the amplitude and the green channel the phase modulation. 
The SLM controller extracted these information and applied the two channels to the respec- 
tive panel. We have performed calibration measurements to ensure that both the amplitude and 
phase modulation exhibited a linear dependence on the applied 8 -bit color value between and 
255. A CCD camera allowed us to record images of laser fields in the Fourier plane of lens 5. 

To conform this experimental section to the conventions introduced in section [3Jl we remark 
that we shaped a set of test fields both in amplitude and phase in the initial plane at z = Zi 
which coincided with the two SLM panels and subsequently minimized the size of a focal 
spot in the target plane at z = Z2 which coincided with the CCD camera chip. For our proof- 
of-principle experiments we ignored polarization effects and considered a set of scalar fields 
Ei{x^y^z) = A/^^^'(^'^'^) (/ = I...N) where the CCD camera detected the associated intensity 
Ii{x,y,Z2) "^ \Ei{x,y,Z2)\^- Both the QMEiIO and the QME:SSO were assembled from these 
fields according to the scalar expressions ([T2l) and (fTTl) , respectively. 

The amplitudes Ai{x,y^Z2) were determined by simply recording an associated intensity im- 
age Ii{x^y^Z2) and subsequently taking the square root. We used the three-step phase retrieval 
algorithm described in Ref. ll2T1l to retrieve the phase modulations 0;(x,}^,Z2). This algorithm 
is based on interference with a reference field £^r(x,j,z) = A^{x^y^z)e^^^^'^'^\ The reference 
field's intensity was distributed uniformly over the ROIs considered in our experiments by 
adding a square phase 0r oc (x^ -\-y^) to the reference field using the SLM phase panel. More- 
over, a constant phase term 0/. = 2nk/3 was added for ^ = 0, 1, 2, and the superimposed fields 
Ei{x^y^Zi ) +£'r(x, j,zi )e^^'' were then encoded onto the SLM. The associated intensity distribu- 
tions explicitly read 



/,• ;^ = |A,f + |ArP + 2 V^4;a^ cos ((/)R - 0, + 2;r V3) 

=: Ibg,i + 7i cos( A0,- +27rV3), (20) 

where the spatial coordinates (x,j,z) were omitted for brevity. These three intensity distribu- 
tions represent a 3-dimensional equation system with three unknowns /bg,/, Jt and A0/. The latter 
is explicitly obtained as 

A0/ mod 2n = atan2{(l - cos(2;rV3)) (//,o -//,2) , sin(2;rV3) (2/,-,i -/,-,o -^^2)} , (21) 

where atan2{^, (^} is the two argument arctangent function corresponding to the argument of 
the complex number <^ +iC- Note that the reference phase 0/? cancels out when calculating 
the operator elements due to multiplication of a complex conjugate field oc ^-i<fe(^))^'^) with a 
complex field oc ^i</>r(^'>^'^) . Therefore Eq. ([2T1) yields the adequate phase modulation required to 
assemble the QME operators. 

During the course of our experiments we verified the linearity of our optical system by per- 
forming a comparison between what we term the "experimental superposition (Exp-S)" and 
the "numerical superposition (Num-S)". The Exp-S refers to the case where the set of QME 
optimized superposition coefficients at is used to encode the optimized superimposed field 
onto the SLM. The CCD camera then detected the intensity lExp-s{^^y^Z2) corresponding to 
this encoded optimized field. The Num-S utilizes the fields Ei{x,y,Z2), which were individ- 
ually measured to assemble the QME operators, in order to numerically determine the in- 
tensity distribution as lNum-s{^jyjZ2) °^ \Y4Li^iEi{x,y,Z2)\ • Crucially, linearity is verified if 



^Exp-s(-^7 J?^2) = ^Num-s(-^7 ^7^2)- This is indeed observed in our experiments as demonstrated in 
the following subsection which features a comparison of experimental and numerical intensity 
distributions. 

4.2. Results and discussion 
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Fig. 10. SLM encoded field modulations and resulting beam profiles, (a) Ring mask RGB 
image as encoded onto the dual panel SLM. (b) Associated Bessel beam created in the 
CCD camera plane, (c) Aperture RGB image as encoded onto the dual panel SLM. (d) 
Associated Airy disk as detected by the CCD camera. The yellow bar in (b) represents 2 
times the SOIM w^ of the Bessel beam's central core, w in (d) is the SOIM of the Airy 
disk. 

In our experiments, we used N =\\ non overlapping amplitude ring masks with a constant 
phase modulation as fields of interest Ei{x^y^z\). After propagation through the Fourier lens 5 
(see Fig.O the resulting fields Ei{x^y^Z2) form a set of Bessel beams. Figure [TOt a) shows the 
largest ring modulation encoded onto the SLM with the resulting Bessel beam shown in Fig. (b). 
As this particular Bessel beam comes along with the highest NA compared to the Bessel beams 
created with smaller ring modulations, the beam shown in Fig. (b) exhibits the smallest central 
spot of all beams realized in our experiments. The SOIM of the Bessel beam featuring the 
smallest core is denoted as wb and used as reference for the measurements presented below. 
For comparison Figure[TOtc) depicts a circular aperture which is encoded onto the SLM in order 
to observe the Airy disk (see Fig. (d)). The SOIM of the Airy disk is approximately L5 times 
larger than the core of the reference Bessel beam as expected [22|. 

The results of the performed QME spot size minimization are shown in Fig. [TT] for dif- 
ferent sizes of the ROI. To begin with, the comparison of the Num-S intensity distribution 
^Num-s(-^7j7^2) (top row) and the Exp-S intensity distributions /exp-s (•^73^7^2) (bottom row) 
clearly reveals good agreement and thus verifies the linearity of our optical system as eluci- 
dated above. For completeness, the central row shows the Exp-S superposition in RGB format 
as encoded onto the SLM. The color code features a blue channel representing the amplitude 
modulation from (black) to 1 (blue) and a green channel corresponding to phase modula- 
tions from (black) to 2;r (green). Next, we conclude from the measured relative SOIM w/wb 
that the spot size decreases if the ROI size is reduced. The reduced spot size is achieved at 
the expense of the spot intensity which is redistributed to a ring outside of the ROI similar 
to the theoretical results presented in section [34l and Fig. [71 Referring to the Exp-S data, for 
R = 7 pixel the spot size is reduced to 72 % of the size of the reference Bessel beam's core and 
even further to 50 % for 7? = 4 pixel. The latter result is somewhat vague, though, due to the low 
spot intensity which may be truncated by the sensitivity threshold of the CCD detector and thus 
may appear smaller. However, our experimental results overall clearly verify the QME concept 
applied to spot size minimization. Moreover, the results strongly suggest that the QME opti- 
mization may indeed squeeze spots to the subdiffractive regime since the optimal superposition 
of Bessel beams not only beats the Airy disk but also the reference Bessel beam. 

We have performed an additional comparison of the experimental results to the theoretically 




Fig. 11. Experimental QME spot size minimization. Top row. /Num-s(-^ 5)^5^2) for different 
ROI radii in pixel as indicated in the top left corner of all graphs shown. The ROI is exem- 
plary indicated as a dashed ring in the left hand side intensity distribution. The number in 
the bottom left corner represents the SOIM w in units of the reference SOIM wb- Central 
row. Optimized experimental distribution as RGB encoded onto the SLM. Bottom row. In- 
tensity distributions /Exp-s(-^5 ^5^2)- The relative SOIM w/wb is indicated in the lower left 
corner. 



predicted ones: The relative SOIM w/wb was evaluated for the Num-S for ROI radii ranging 
from R=\ pixel to 7? = SOpixel in steps of one pixel (see Fig. [12]). In comparison to the cor- 
responding graph shown in Fig. [8] we observe agreement in terms of both the general decrease 
of w/wb with decreasing R and the peaks occurring periodically along the /?-axis. Surprisingly, 
for R> 10, we do not observe values of w/wb ^ 1 ; that is the QME spot minimization does not 
deliver, as one might naively expect, the reference Bessel beam featuring the smallest SOIM 
wb of all Bessel beams considered. This is due to the ring structure of Bessel beams which 
significantly adds to the SOIM and therefore to the spot size in the case of large ROIs. As a 
consequence, the QME optimization aims to superimpose the set of Bessel beams in a man- 
ner that the rings are destructively interfered within the ROI, only retaining the central core. 
Given that the ring structure is essential for the reduced core size of Bessel beams compared to 
the diffraction limited Airy disk we observe values of w which are larger than wb . Indeed, the 
values of w are fairly close to the size of the Airy disk which was determined as w ^ 1.5wb. 
This strongly suggests that for large ROIs the QME spot minimization yields a superposition 
of Bessel beams which closely resembles the diffraction limited Airy disk. Finally, we remark 
that for R<1 the relative SOIM w/wb becomes very small and has to be taken with care due 
to the possible truncation of the spot mediated by the CCD detector's sensitivity threshold as 
mentioned above. 
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Fig. 12. Relative SOIM w/w^ within the ROI depending on the ROI radius R in pixel for 
numerical superposition of the fields Ei{x,y,Z2) in the CCD plane. 



5. Discussion and Conclusion 

We have theoretically and experimentally demonstrated a novel approach based on quadratic 
measures eigenmodes that enables the optimization of different optical measures. The theory 
that we employ is rigorous and based on considering the light-matter interaction as a quadratic 
measure originating from the fields described by Maxwell's equations. Excitingly, we can de- 
fine many quadratic measure operators to which our approach is applicable (see Table [T]). The 
method is thus very powerful and the generic nature of our approach means that it may be ap- 
plied for example to optimize the size and contrast of optical dark vortices, the Raman scatte- 
ring or fluorescence of any samples, the optical dipole force, and the angular/linear momentum 
transfer in optical manipulation. In the present paper we have verified the rigor of the method 
by demonstrating experimental spot size operator and intensity operator optimization using 
Laguerre-Gaussian and Bessel light modes using a dual SLM approach to implement the tech- 
nique. We envisage the QME approach as providing a powerful and versatile theoretical and 
practical toolbox. Our generic approach is applicable to all linear physical phenomena where 
generalized fields interfere to give rise to quadratic measures. 
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